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Chapter 1 


INTRODUCTION 


1.1 Introduction to the Thesis 

With an ever-rising demand for superior quality of steel at a competitive price, steel in- 
dustries globally are on the look out for processes with the following features for a rapid 
attainment of end chemistry: 

(i) Faster alloy dissolution to reduce processing time, 

(ii) Higher alloy recoveries to have an edge on cost and cleanness front and 

(iii) minimal thermal and particulate concentration gradient in the melt. 

It can therefore be said that mass transfer and liquid mixing phenomena are perhaps the 
most critical operations in terms of quality steel production. 

The melting and/or dissolution of solid additions in liquid steel melt is an integral 
unit operation in steel making industries. Large tonnage of ferro alloys and deoxidizer 
elements are added to liquid steel bath on a routine basis to adjust the chemistry of the 
melt and thus, produce a grade of steel to meet the customer requirements. Melting 
and dissolution of additions consume a large amount of energy and have an important 
bearing on the overall process economies. Quantitative frame works are therefore needed 
so that these can be predicted and also important informations on energy consumption, 
complete melting/dissolution times, optimal size ranges vis-a-vis the processing time etc. 
are known a-priori. This information, can help operators to run the melt shop more 
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economically and effectively. 

In order to ensure a sound solidified product for customer, it is important that virtu- 
ally no stratification of heat and mass takes place in the melt prior to casting. There has 
been a constant endeavour to exacerbate the rates of these transport process with obvious 
repercussions in the steel industries and thus, justifiably, these have been the subject of 
scientific research for more than two decades. Homogenization of melt with respect to 
composition and/or temperature is the fundamental requirement to guarantee a sound 
casting. Homogenization of steel prior to casting may be attained most economically by 
stirring liquid steel in a ladle (a cylindrical vessel with a slight taper towards bottom) in 
several ways. The most common configuration is bottom gas injection through a centre, 
or eccentric porous plug, or lance or tuyers. Liquid mixing in a gas-stirred bath occurs 
by the recirculatory and turbulent motion of the fluid induced by the injected gas. It 
is therefore important to develop a good understanding of the induced turbulent fluid 
motion in conjunction with the dispersion or homogenization process. 

Mixing time is typically defined as the time within which, a vessel’s liquid can reach a 
state of chemical and/ or thermal stability. Mixing phenomena should investigated both 
theoretically and experimentally, but, owing to the large size of industrial processing 
units, high operating temperatures and visual opacity of liquid steel, it is practically im- 
possible to. carry out any meaningful experimental programme. Mathematical modelling 
provides a reasonable alternative to study the mixing phenomena in industrial processing 
units. 


1.2 Previous work 

1.2.1 Liquid mixing phenomena 

So far as mathematical modelling of mixing times is concerned, different approached have 
been adopted by researchers. A summary of various procedures available for calculations 
has recently been presented by Mietz and Oeters [1]. From a conceptual stand point, two 
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different approaches have been applied and these include, modelling via numerical solu- 
tions of the species conservation equation and the circulation time model. In conjunction 
with mathematical modelling, experimental measurements of mixing time in water mod- 
els and industrial argon stirred ladles were also carried out. Such measurements were 
often applied to validate model predictions. 

In presence of a two dimensional velocity field, the mass conservation of an inert tracer 
i (e.g., TTii) can be expressed (in a cylindrical polar co-ordinate system) via the following 
convection-turbulent diffusion equation: 


S 6 . .16, . 5 , 5mi. 1 6 , Snu. 

f 

The eddy diffusivity Dj^ eddy kinematic viscosity, (=^), are conventionally 

taken to be numerically equal. It is therefore apparent, that provided the flow parameters 
(u and v) and turbulence quantities {he) are known with reasonable certainty, mixing 
time [viz., mi{z, r, t) fields] can be fairly accurately predicted because physically realistic 
and accurate boundary conditions can be applied to Eq. (1.1). This has been a popular 
approach in numerous mathematical modelling studies in which mixing phenomena have 
been studied in aqueous as well as industrial gas stirred ladle system. 

In a much earlier study, Szekely et al. [2] reported a combined theoretical and experi- 
mental investigation of liquid mixing phenomena in industrial scale argon stirred ladles. 
Mixing rates in three different size industrial ladles (viz., 7, 40, and 60 t) were measured 
for various operating conditions by monitoring a transient local concentration of a ra- 
dio active tracer addition. In conjunction with their experimental work, they predicted 
mixing times through solution of Eq. (1.1) via a two dimensional elliptic, turbulent flow 
model. Despite some simplifying assumptions, experimental measurements and numeri- 


cal predictions were found to be in a reasonably good agreement. Similar observations on 
water models of argon-stirred ladle have also been reported by Ying et aJ. [3], Salcudean 
and co-workers [4], Mazumdar and Guthrie [5] and more recently by Joo and Guthrie 
[6]. Computational studies of mbdng phenomena in industrial size and water model la- 
dles confirm that gas flow rate and vessel geometry affect mixing considerably and that 
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mixing rate vary significantly from one location to another in such gas stirred vessels. 

In subsequent investigation El-Kaddah and Szekely [7] improved their earlier hydrody- 
namic model to predict the turbulent flow fields in 6 and 40t industrial ladles. Embody- 
ing the predicted flow parameters and turbulence quantities in the species conservation 
equation and invoking some thermodynamic equilibrium relationship, desulphurisation 
kinetics in the industrial argon stirred ladle were numerically predicted and reasonably 
good agreement with experimental measurements was demonstrated. 

Mixing phenomena in gas stirred ladles designed for single and duel plug bubbling op- 
erations was reported recently by Joo and Guthrie [6]. Their computational results and 
experimental observation confirmed that, as the bubbles are moved off-centre, angular 
momenta increase, reducing mixing time considerably. It was found that a mid-radius 
placement of a porous plug represents an optimum location for single plug bubbling, 
while diametrically opposed, mid radius placement of bubblers were found to be best for 
dual plug bubbling. These numerical and experimental modelling exercise indicate that 
mixing phenomena in gas stirred ladle systems can be reasonably accurately predicted 
from first principles. 

An alternative to the differential modelling approach to investigate mixing phenom- 
ena in gas stirred ladle systems was first applied by Sano and Mori [8]. The calculation 
procedure, commonly known as the circulation time model, is essentially based on the 
concept that the circulation time is proportional to the mixing time. Through an energy 
balance (e.g., at steady state, the rate of potential energy supplied by the rising bubbles 
is balanced by turbulence energy and interphase frictional dissipation losses), Sano and 
Mori derived an expression for estimating circulation time in terms of operating variables. 
Furthermore, by assuming mixing time to be equivalent to three times the circulation 
time (e.g., equivalent to a degree of mixing equal to 100 ±5%), mixing time in gas stirred 
cylindrical vessels containing molten iron were predicted mathematically and compared 
with experimental measurements. The calculated results were however, found to agree 
only roughly with measurements. In a later study, the same approach was also adopted 
by Stapurewicz and Themelis [9] to investigate mixing phenomena theoretically. 

A similar approach was adopted by Krishnamurthy et al. [10] to theoretically inves- 
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tigate mixing phenomena in gas stirred baths. On the the basis of many experimental 
data derived from water model studies, Krishnarnurthy et al. suggested that circulation 
number (= circulation time/ mixing time ), in contrast to the earlier supposition of Sano 
and Mori [8], is not a constant and can rather assume a value between 2 and 12 ( for any 
given degree of homogeneity). In a subsequent study, Krishnarnurthy [11] applied a ’’tank 
in series model” to compute the circulation number and hence, the mixing time. Model 
prediction were assessed against experimental measurements and reasonable agreement 
claimed. However, for high temperature industrial systems agreement was found to be 
less satisfactory. 

The models of mixing time adopted by Sano and Mori [8], Stapurewicz and Themelis 
[9] as well as Krishnarnurthy and coworkeis [10,11], are some what simplistic. Although, 
through tuning of some key parameters in the model equations, the circulation time 
models can be made to describe a particular set of experimental conditions, nevertheless, 
such empirical models, in general will have only limited utility, since these cannot be 
readily extrapolated with much confidence to other configurations. The concept applied 
in the formulation of the model equations are fundamentally weak since these do not 
incorporate the physics of the mixing process (e.g., that macro-mixing process via bulk 
convection, turbulent diffusion and molecular diffusion phenomena ) adequately. Thus, 
Mietz and Oeters [1] suggested that the recirculating flow, turbulent diffusion from the 
tracer enriched cloud and the mass exchange between the dead zone and the remaining 
volume of the liquid in tanks are parameters that have to be considered in such ’’tank in 
series” or ’’recirculation time” models. 

1.2.2 Mass transfer between solid and liquid 

Thermal and mass interactions between solids and fluids are characteristics of many pro- 
cesses (viz., alloying addition, powder injection and so on ) that are carried out in refining 
ladles. Over the years, several experimental studies have been reported and new empiri- 
cal heat and mass transfer correlations have been developed to estimate melting and/ or 
dissolution rates in gas stirred ladles. These investigations, in general, have appeared to 
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indicate that in turbulent metal processing systems, given that the intensity of turbulence 
is often appreciable, classical correlations for translatory flow are insufficient to connect 
the hydrodynamic transport phenomena to process rates. Thus, in order to estimate 
heat and mass transfer rates in gas stirred systems effectively, functional relationships 
embodying the combined influence of flow and turbulence have been proposed. 

Szekely et al. [2] investigated the the dissolution rates of graphite rods dipped into 
industrial size gas stirred ladles. To estimate the mass transfer coefficient theoretically, 
the correlation proposed by Lavender et al. [12] was applied e.g., 

Sh = 2 + Q.72Re'^-^Rel^Sc°-^ (1.2) 

A comprehensive theoretical and experimental investigation on melting phenomena 
in turbulent recirculatory gas bubble driven aqueous and high temperature systems has 
been reported by Taniguchi and co-workers [13,14]. These authors measured the melting 
rates of ice sphere using a photographic technique and compared them with estimates 
derived from three diff'erent heat transfer correlations reported in literature. Such com- 
parisons demonstrated explicitly that the classical Ranz-Marshall correlation ( viz., 
Nu — f +0.6(i?e)°'®(Fr)°'^^) underpredicts heat transfer rates in gas bubble driven sys- 
tems by about 40% to 50%. On the other hand, the heat transfer correlation proposed by 
Whitslker [15] (e.g., a modified version of the Ranz-Marshall correlation and applicable 
to turbulent flow situation) viz., 

(Nu - 2) = [0.4Re° ® + 0.06Re*’ ®®]P^° '‘ ( -^) ” (1-3) 

\ muo/ 

was found to produce better estimates and was therefore more appropriate and adequate 
for investigating melting phenomena in both aqueous and high temperature gas stirred 
baths. In a later study, Szekely, Gravet and El-Kaddali [16], proposed a correlation 
(applicable in the range of values for Re (between 100 and 2000) and Tu (> 0.15) ) to 
estimate melting rates of ice cylinders in gas bubble driven system according to: 

Nu = 0.338(Pe ru)‘’V°'^^ 


(1.4) 




1.2 Previous work 


7 


111 order to estimate various non-dimensional numbers and hence the corresponding melt- 
ing rates, a photographic technique in conjunction with a turbulent flow model was ap- 
plied. 

A similar approach, was adopted by Mazuradar and coworkers [17,18], for estimating 
mass transfer coefficients for the dissolution of benzoic acid cylinders in a gas stirred 
aqueous bath. On the basis of experimental measurements and numerical computations, 
it was demonstrated that mass transfer coefficients for vertically submerged cylinders in 
a water model ladle can be adequately expressed via a similar correlation of the type 
expressed via Eq. (1.2), e.g.. 


Sh = Q.73Re^-^Ref^^Sc'^-^^ (1.5) 

in which the nominal and turbulent Reynolds numbers {Reioc,r cmdRet) were the local 
flow variables and defined as ( ^ -- ^- --) and (^), respectively.' In those studies, the 
mean and fluctuating velocity scales embodied in the definition of nominal and turbu- 
lence Reynolds number were estimated theoretically from the k — e model of turbulence 
(e.g., «=\/0.66 k). In a separate study, Mazumdar et al. [19] extrapolated the above 
equation and demonstrated that isothermal dissolution rates of steel cylinders in a car- 
bon saturated iron melt can be predicted with reasonable certainty via Eq. (1.5). These 
investigations, together with many others, are therefore indicative of the fact that appro- 
priately designed water model investigations can be reasonable and effective alternatives 
for studying various high temperature phenomena in ladle refining operations. 

Solid and particulate dissolution rates in laboratory scale, gas stirred iron-carbon 
melts were reported by Wright and coworkers [20,21]. Mass transfer rates for the isother- 
mal dissolution of steel rods were measured in two different sized (1 and 25kg) baths for a 
wide range of gas flow rates. It was confirmed experimentally that gas injection enhances 
the rate of dissolution considerably. In their subsequent study, recarburisation phenom- 
ena in a gas stirred bath were investigated via multi-particulate dissolution experiments. 
To compare measured and predicted particle mass transfer rates, Wright and coworkers 
applied the following correlation, derived from Kolmogoroff’s theory of local isotropy, e.g., 
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\ 0.25 

( 1 . 6 ) 

It was demonstrated that Eq. (1.7) can adequately simulate multi particle dissolution 
phenomena (e.g., recarburisation etc.) in high temperature gas stirred system. 

Asai and coworkers [22] have compiled and discussed some of the reported mass 
transfer correlations from the viewpoint of their applicability to ladle refining operations. 
One of these correlations, viz., 

Sh = 0.079 Re°-^ (1.7) 

was subsequently applied by Koria [23] to estimate liquid velocity (e.g., Re) from the 
measured dissolution rates (e.g., Sh) of solid oxalic acid compacts. The fluid velocities 
estimated from the measured dissolution rates were compared with corresponding theo- 
retical estimates obtained from simplified macroscopic relationships and found to be in 
excellent agreement. This appears to contradict the findings of many other investiga- 
tions, since object Reynolds number alone has been shown to be inadequate to simulate 
observed heat and mass transfer rates in gas stirred ladles. Nevertheless, in view of the 
relatively small size of the vessel (L=0.15m and D= 0.15m )applied in conjunction with 
low gas flow rates, it is not unlikely that intensity of turbulence in that system was small 
and therefore, a correlation (essentially valid under laminar flow conditions ) embodying 
only object Reynolds number produce reasonable estimates of dissolution rates. It is also 
useful to note that the intensity of turbulence within a gas stirred ladle is strongly spa- 
tially dependent. Furthermore, flows across the ladle’s base and regions deep within the 
ladle are among the slowest. Consequently, the regions in which dissolution is measured 
must greatly affect the conclusions drawn. 

The subsurface motions and trajectories of solid additions to gas stirred ladle systems 
were investigated both theoretically and experimentally by Mazumdar and Guthrie [24]. 
Their study suggested that under industrial conditions, buoyant additions, irrespective of 
their size, have practically no chance to undergo prolonged subsurface motion but rather 
almost instantaneously resurface. Similarly, they found that heavier additions would al- 
ways sink to the bottom of the vessel. In contrast to these, neutrally buoyant additions 
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(e.g., some grade of ferromanganese), have the potential to undergo prolonged subsurface 
motion, and therefore melt within the bulk metal bath. On the basis of their physical 
and mathematical modelling studies, it was suggested that buoyant additions such as Al, 
Fe-Si, etc. would, in conjunction of variable amounts of carry-over slag, exhibit poor and 
erratic recovery rates, if added to ladles during argon bubbling operations. In their con- 
text, these authors highlighted the role of the C.A.S (Composition Adjustment through 
Sealed argon bubbling) procedure. A summary of various correlations related to heat 
and mass transfer have been presented in Table (1.1) 


1.3 Scope of Present Work 

Two phase modelling of mixing phenomena in materials processing has few and far be- 
tween. In particular, combined Lagrangiaii-Eulerian modelling approacli has not been 
applied so far for predicting mixing rate in gas-stirred reactors. The predictions from 
class of models in general depend on empirical inputs (such as, bubble diameter, drag 
coefficient vs Reynolds number relationship etc.) and since these are somewhat uncer- 
tain, particularly in high temperature gas stirred baths, it is important to assess the 
influence of these empirical parameters on the rate of predicted liquid mixing. In com- 
bined Lagrangian-Eulerian approach, the motion of gas phase is computed numerically 
in a Lagrangian frame of reference in the form of steady stream of bubbles while the 
liquid motion together with fluid turbulence is calculated by Eulerian scheme [26]. The 
steady state, two phase, fluid flow model in conjunction with the standard /c-c model of 
turbulence is used in the present work for simulating the flow phenomena. These models 
are coupled with an unsteady state, species conservation equation for predicting mixing 
times. For mathematical analysis of mixing phenomena, it is important that predicted 
results are independent of grid size distribution and incremental time step. So, in the 
present work influence of grid size and time step on predicted mixing times has been 
investigated computationally. 
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Table 1.1: A summary of various correlations applied to investigate heat and mass trans- 
fer phenomena in metal processing operations 


Correlation 


Investigator 


Remarks 


Sh = O.rSRe^-^^Rcf^^Sc^-^^ 


Mazumdar, 
Khajani and 
Ghosh [17] 


developed for 
analysis of 

dissolution phenomena 
in gas stirred bath 


Sh = 2 + 0.7Re“®ReO-255c‘’-33 


Szekely, 
Wang and 
Lehner [2] 


for spherical 
geometry in 
turbulent shear flows 


{Nu - 2) = [0.4ii!e° ® + 0.06Re° ®®] 




Taniguchi 
et al. [13,14] 


for spherical 
geometry 


Nu = 0.8(Reru)°-8Pr°-33 


Szekely, 
Gravet and 
El-Kaddah 
[16] 


for analysis of 
melting phenomena 
in turbulent gas 
stirred reactors 


Sh = 0.07ReO'^5'c°-356 


Koria 


for analysis of 
dissolution phenomena in 
ladle metallurgy 
steelmaking operations 


Sh = 0.323/t:(;°®Pc'’-'’^ 


Geiger and 

Poirier 

[25] 


Forced convection 
over a flat plate 
valid for most fluid 
including liquid metals 
(Approx. Integral Technique) 


Sh = 0.646Re°®5c°-“ 


do 


The average mass-transfer 
from a flat plate to a fluid 
(Approx. Integral Technique) 


Sh = 0.332Re“®5c°-3'‘5 


do 


The local mass-transfer for 
transfer from a flat plate 
to a fluid is 

given by (Exact Solution) 


Sh = 0.664Pe°®5c‘’-^''^ 


do 


The average mass-transfer for 
transfer from a flat 
plate to a fluid is 
given by (Exact Solution) 
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Significant vaxiations in local mixing rates in different regions of the gas-stirred vessel 
were reported by some investigators. In contrast some researchers have reported con 
tradictory observations. Thus, in the present work, the position dependency of mixing 
in gas-stirred vessel has been analysed rigorously through mathematical model studies 
for different size vessels to identify the reasons for the contradictory observations. The 
reliability of the present model has been evaluated by comparing model predictions with 
experimental data. In this regard, the variation of mixing times with liquid depth has 
been investigated computationally. 

Several mass transfer correlation have been reported (see section 2.2) in the literature 
which are essentially developed by curve fitting technique. These correlations are nor 
rnally expressed in dimensionless form by using characteristic dimensionless groups viz-, 
Sh, Re (local), Re (turbulent), and Sc numbers. Since, for forced convective mass trans 
fer, Sherwood number (^) can be expressed as a function of local Reynolds number 
and turbulent Reynolds number (^), an attempt has been made in the 
present work to train these correlations by an artificial intelligence teclinique namely 
Genetic Adaptive Search (GAS), using the experimental and theoretical data reported 
in the literature so as to deduce optimal values of pre-exponent and exponents in a 
particular mass transfer correlation. Finally, a steady state two phase turbulent flow 
calculation procedure has been applied to predict mass transfer rates from first principle 
in two different gas stirred system and assess the adequacy and appropriateness of the 
GAs optimized models. 


1.4 The Layout of the Thesis 

The main body of the Thesis consists of five chapters. Chapter-2 deals with numerical 
predictions of mixing time in a gas-stirred bath. In this chapter, a steady state two 
pluise turbulent flow model and an unsteady state, species conservation equation ( t 
were solved for predicting mixing time) are explained. In Chapter-3, the developmen 
mass transfer correlations by applying Genetic Algorithms has been presented. In 



1.4 The Layout of the Thesis 


12 


chapter modelling of mass transfer from the first principle has also been briefly discussed. 
Cha.pter-4 summarises the general conclusions derived from the present work. Finally, 
Chapter-5 concludes this Thesis with the recommendations for future work. 



Chapter 2 


NUMERICAL PREDICTION OF 
MIXING TIME IN A 
GAS-STIRRED BATH 

2.1 Assumptions in Modelling 

A combined Lagrangian-Eulerian fluid flow model in conjunction with appropriate species 
conservation equation has been applied in the present work to theoretically predict mixing 
phenomena in gas-stirred vessels. Several assumptions and idealisation are summarised 
below: 

(i) The presence of any overlying second phase is ignored and melt-air surface was as- 
sumed to be flat. 

(ii) The melt is assumed to be isothermal. 

(iii) Following axial symmetry in the present geis injection configuration, variations of 
flow properties along the angular direction are ignored. 

(iv) Discrete mono-size bubbles are assumed to be generated at the lance tuyer/tip. 
Thus, the size of the bubbles created at the tip is known a-priori. The bubble frequency 
is calculated from the principle of volume continuity. 

(v) Bubble - bubble interactions are assumed to be negligible and hence drag coefficient 
correlation, applicable to the single bubble-fluid system is applied to the appropriate 
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range of Etvos and Reynolds number. 

(vi) Bubbles within the upwelling plume contribute to the generation of turbulence and 
this in turn affects the bubble motion. 

(vii) The fluid is essentially Newtonian and incompressible and finally, 

(viii) The flow is assumed to be steady, turbulent and recirculatory. 


2.2 The Mathematical Model 


2.2.1 Governing equation of liquid mixing 


In the presence of a two - dimensional velocity field, with no generation, the conservation 
of mass is expressed in terms of cylindrical polar co-ordinates as 


+ —(aipurm,) = (2-1) 

Where r, effective exchange coefiicient (//<) can either be approximated using an effective 
viscosity formula proposed by Sahai and Guthrie[26], or from an appropriate turbulence 
model. 


2.2.2 Governing equation of liquid phase motion 

The knowledge of flow fields (i.e u and v ) and turbulence parameter (i.e // ), are required 
to solve the concentration (mass) equation. 

A brief outline of the combined Eulerian - Lagrangian two phase, two dimensional tur- 
bulent flow model developed earlier by Mazumdar and Guthrie[26] is presented below. 
The governing equations together with their boundary conditions applicable to axisym- 
metrical gas stirred ladles are written in terms of cylindrical co-ordinates. , 

Equation of continuity 


2.2 The Mathematical Model 


5 1 j 

Y^{piaiu) + --[p^aiTv) = 0 


Equation of motion in the axial direction 


l^(a,p,uu) + llipmruv) = -a,| + + 

Equation of motion in the radial direction 


^{(Xipiuv) + ]^^{pioiirvv) = 


(5p . 5 , 6v. 6 , Sv. 




2.2.3 Governing equation of liquid phase turbulence 

In the present work a modified form of standard coefficient k — e model of Launder and 
Spalding [27] was used. The governing equations for turbulence kinetic energy, k and its 
dissipation rate, e, can be represented in the cylindrical polar co-ordinates as 
Equation of turbulence kinetic energy, k 


f^(pmpK) + ll(p,mrvK} = ^ 


+o:j(^fc “ Pl^ + 


Equation of dissipation rate of turbulence kinetic energy, e 


Pt Se. 


PrSe. 


ip.a,pe) + + 

,CiGtl CjFie" , 

Ctl{ + ) 

K K K 
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In Eqs(2.5) and (2.6), Gk, the volumetric rate of turbulence production, is defined as 


Gk = Pr 


/(5uV f Sv\^ / 


'"(!« Sv\^ 
Sr ^ Sz) 


The eddy viscosity, He, appearing in Eqs (2.3) and (2.4) is defined as 


f^e — IJ^ + fJ>L 


Where, 


(2.7) 


( 2 . 8 ) 


/ir = 


C^pik'^ 


(2.9) 


The five empirical constants of the /c — e model viz., Ci,C 2 , Cfi,ak, and <t, were assigned 
standard values as summarised in Table (2.1): 


Table 2.1: Values of constants in the « — e turbulence model 


Constants 

Cl 

C2 

(^K 


Cd I 

Values 

1.44 

1.92 

1.0 

1.3 

1.0 1 


2.2.4 The equation of motion of the gas phase 


The analysis of motion of the gas phase has been carried out via a Lagrangian approach, 
in which the trajectories of individual bubbles are stochastically determined in space by 
solving two simultaneous ordinary differential equation in time. The Newton’s second law 
of motion can be represented along the two co-ordinate axis, by the following equations : 
Along axial direction z 


(1 + 0 . 6 ^)^ = -U) + 1 . 5 ^ 

pg^ dt Uipg Pg 


oz or 


+ (l--f)s (2.10) 
Pa 
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Along the horizontal / radial direction r : 


(1+0.5^)^!^ 
Pg at 


3piCoRe 


(Vb, 


V) + 1.5^i 


ut + V^'’ 


( 2 . 11 ) 


The two corresponding kinematic relationships which define the trajectories of bubble’s 
within the two dimensional flow domain are : 


% = Vb. ( 2 . 12 ) 

and, 

The instantaneous drag coefficient, Co, in the Eqs. (2.10) and (2.11) is evaluated by us- 
ing an experimentally determined correlation for oblate spheroid / spherical cap bubbles, 
and expressed as 


Co = 


0.622 
-f- 0.235 


for-, (500 <Re< 5000) 


(2.14) 


Furthermore, Re in Eqs. (2.10) and (2.11) is the Reynolds number, and is estimated on 
the basis of relative velocity vector Ur, i.e 


Re = pjdft 


Pb-Ufl 


t^i 


(2.15) 


In eqs. (2.10) to (2.13), Vb,z and VB,r are the instantaneous components of bubble motion 
in the axial and the radial directions respectively. U and V, the instantaneous component 
of fluid motion, are estimated by summing up the relevant mean velocity components (u 
and v) with the corresponding fluctuating velocity components, assuming that turbulent 
fluctuations are essentially isotropic and follow a Gaussian probability distribution func- 
tion. For example, 


U = n + ,l>(jy ( 2 . 16 ) 

Where, V’ is a normally distributed random variable and varies between -1 and 1. Thus, 
during the rise of bubble through the liquid, the velocity fluctuations are estimated from 
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(i) the predicted k field and 

(ii) the random variable , 0 

I he velocity fluctuations are assumed to remain constant for a time span, equivalent to 
the eddy life time . Furthermore, the life time of the eddy, during which a characteristic 
velocity fluctuation is assumed to exist, is estimated from the turbulence theory accord- 
ing to 


r = 0.195- (2.17) 

e 

For a scalar control volume the corresponding gas voidages can be estimated from the 
principle of volume continuity according to ; 

a = volume of gas ^ y ^ ^2 is) 

^ volume of control volume AKw N V • ; 

Where, ^ represents the mean residence time of N bubbles in the control volume. 

The sum of the residence time terms can be expressed as: 


( 2 , 19 ) 

Hence, the gas volume fraction ag, for a particular trajectory and specific to a control 
volume, can be represented as: 


^ Qt{N + l)At ^ 20) 

" 2AVc,„ ^ ^ 

Once ag values have been determined for each of the control volumes, corresponding 
liquid volume fraction oj, can be easily determined, from -I- aj = 1.0. 

It is now well accepted that the bubbles in the two phase plume region contribute sig- 
nificantly to the production of turbulence. The additional turbulence produced by the 
bubbles can be considered to be equivalent to the shear work performed on the liquid by 
the bubbles and therefore, the corresponding volumetric flow rate of turbulence genera- 
tion , Pb, can be represented as: 
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P — n YJ \2jj. 


( 2 . 21 ) 


Where, Vg is the resultant bubble velocity and Uf is the resultant fluid velocity . Cf, in 
eq. (2.21) is a constant and takes a value between 0 and 1. The value of Ct depends on 
the bubble population, bubble size etc. in the two phase plume region. The constant Cb 
is normally determined by trial and error. 

The momentum exchange terms and Fr can be evaluated by assuming that the drag 
force experienced by bubbles acts with equal magnitude but in opposite direction. The 
expression for F^ and F^ can be represented in the compact tensorial notation as: 


Fi = 


Qi 


NAV„ 


m=n 

SJ. 


*r,m SfilCgRe 

4d? 


{VB,i - Ui)dt 


( 2 . 22 ) 


2.2.5 The boundary condition 


The boundary conditions applied to the numerical solution of the set of governing equa- 
tions for species conservation, fluid flow and turbulence in axisymmetrical gas stirred 
ladlfes are as follows: 

(i) At axis of symmetry : 


. = 0; 0<.<L; „ = 0,^ = f=^=0ani^ = 0 (2.23) 


5r 


(ii) At the free surface of the liquid 


z = L\ 0 < r < i?; 
(iii) At the side wall : 


5u Sk Se ,, , Srni 

ti = 0, = — = — = 0 ond — = 0 

6z Sz 6z oz 


(2.24) 


u = v 


= k = e = 0 and = 0 (2.25) 

dz 


r = R] 


0 <z < L] 
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(iv) At the bottom surface : 

z = 0; 0<r<i2; u = v = k = e = 0 and = 0 (2.26) 

or 

The description of boundary conditions is shown in the Fig.( 2.1). 


2.3 Numerical Solution Procedure 

The governing equations have mainly three terms viz., convection term, diffusion term 
and source term. The first order derivatives are convection terms and the second order 
derivatives are diffusion terras. Thus, to numerically solve the set of partial differen- 
tial equations represented by equation (2.1) through (2.13) together with the boundary 
conditions, a control volume based finite difference [28] procedure has been adopted in 
the present work. In this, the flow domain has been discretized into a number of non- 
overlapping control volumes and the governing equations integrated around all these 
control volumes to result into a system of algebraic equations. The staggered control vol- 
umes have been considered for the velocity components (axial and radial velocity ) and 
scalar control volumes for K,e,p and mj. The two-dimensional discretization algebraic 
equation can be written as: 


Oj>$p = ttE^E + O-W^w + O-N^N + + b, (2.27) 

in which, ap is the centre point coefficient while aE,a,w,aN,a,s the neighbour point 
(2 Dimensional) coollicionts, representing the combined influence of fluid convection and 
diffusion and b is the source term. If unsteady term is present in governing equation, 
two terms ap and $p will come into the picture. Then, 

^ P£^y (2.28) 

(2.29) 


b = ScVol. + a°p^°p 
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Figure 2.1: Schematic of the axi-syinrnetric gas-stirred bath illustrating the co-ordinate 
system and the boundary condition used for numerical simulation 




NORTH 



Figure 2.2: Various types of control volumes and grid distributions used in the numerical 
solution scheme 
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Up = oe + Otv + O/V + 05 + a?, - Sp Vol (2.30) 

Here and p°p refer to the known values at time t, while all other values ($p, #Ar, 
and so on ) are the unknown values at time t + At. For steady state governing equa- 
tion the value of Op will become zero. Further, Sc represents the positive component 
of the linearised source term, while the corresponding slope is included in the op term. 
The set of algebraic equations thus derived has been solved via the Tri-diagonal Matrix 
Algorithm (TDMA), through a line by line solution technique. The pressure-velocity in- 
terlinkage in equations (2.2) through (2.6) has been solved by an implicit finite difference 
procedure known as SIMPLE [28] (Semi-Implicit method for pressure linked Equations). 
Since an iterative calculation procedure has been adopted to numerically solve the gov- 
erning partial differential equations, a convergence criteria has been prescribed. To this 
end a residual source has been defined as: 


Residual source = (^ «nfc^’n!i + Sc-ap^p) <10-® {Slunit) (2.31) 


In which, o„6 is the neighbouring point value. The modulus of this residual source has 
been estimated at each nodal point and added up over the entire flow domain for eacli 
individual variable involved. The steps involved for the prediction of mixing time are: 

(i) Velocity field (u, v) and effective viscosity {pe) were calculated from flow and turbu- 


lence models . 

(ii)By using u, v and normalized concentration profiles f(z,r,t) were obtained. Nor- 
malized concentration is expressed as: 

Normalized concentration = S c^en&aS 

Average or bulk concentration was calculated by using the following relationship: 


f„^f^^mi(27rr)drdz 

/o'^//'(27rr)drdz 

In other form, 

- y:'ZmiXVOL{I,J) 

" j:j:voLii,j) 


(2.32) 


(2.33) 
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(iii)Nonnalized concentration profiles (i.e. dimensionless concentration vs. time plots) 
were obtained in the different regions of the vessel. This provides the idea of the rate 
of mixing at various locations in the vessel. The flow circuit of the computer program 
applied to the present study has been iUustrated in Fig. (2.2). 



Figure 2.3: Flow chart of conrputer nrodel used in the present study 
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2.4 Results and Discussion 

2.4.1 The spatial dependency of mixing : the definition of bulk 
mixing time 

The mixing phenomena in a gas stirred bath is related to the turbulent recirculatory fluid 
motion developed due to injection of gas. Since the flow and turbulence are spatially de- 
pendent, it is important to study the variation of mixing rates in different regions of the 
vessel. Normalised concentration profiles have been plotted in Fig. (2.4), as a function of 
time for three different regions ( A, B and C ) of the gas stirred vessel (L=0.80m and R= 
().56m). The rate of mixing and thus mixing times are different for different regions of 
the vessel, 'rims, to quantify a reactor mixing time, a 95% bulk mixing time criterion has 
been applied. This is defined as the time by which the vessels content at every location 
is at least 1 db 0.05 of the equilibrium concentration. Obviously then, 95% bulk mixing 
time is equivalent to the 95% local mixing time of the slowest mixing region/location. 
'Fable (2.2) represents the values of mixing time at three regions of for a axisyminetric 
gas stirred vessel. The region A, which is bottom corner portion of a vessel has maximum 
value of mixing time. In this region species reach at last so this region is considered for 
the evaluation of mixing time. The region B is the central while C is the upper corner 
location of the vessel. It is seen that region B and C have less value of mixing time than 
that of A. As seen, from Table (2.2) and Fig.(2.4), significant variation in local rates of 

Table 2.2; Values of numerically predicted mixing times at different locations of a gas 
stirred bath 


SL. No. 

Location, 

(r,z)cm 

Mbdng time 
(second) 

Region 

region 

1 

(52.5,1.43) 

78 

A 

2 

(52.5,4.29) 

70 

A 

3 

(26.8,38.6) 

59 

B 

4 

(54.8,78.5) 

63 

C 


mixing is observed for different regions 


of vessel. It is found that region A close to the 



Local concentration 
Bulk concentration 
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Figure 2.4: Variation in local mixing time for three different regions, A, B and C in a gas 
stirred vessel ( grid Size = 30 x 26, L = 0.8m, R = 0.56m, dj = 0.015m, gas flow rate = 
6.67 X 10“'* w? Js and Cj = 0-0) 
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bottom of the ladle exhibited the slowest rate of mixing. Hence, it is cleau: from numerical 
predictions that by the time region A approached the 95 percent mixing mark, the bulk 
of the liquid was practically homogeneous. Thus, region A is better representative of the 
bulk 95 percent mixing time of the system than those of locations B and C. 


2.4.2 Influence of grid size distribution on predicted mixing 
time 

The solutions of the governing differential equations for the prediction of mixing time 
are expected to be sensitive to the nodal configuration applied to the numerical solution 
scheme. It is therefore essential to obtain grid independent solutions prior to analysis of 
various aspects on mixing time (viz., influence of time step, location , bubble diameter 
etc). 

To test these, mixing times (via concentration profiles) in a gas stirred cylindrical vessel 
(L=0.80m and 11=0. 5Gm) have been predicted by applying four different grid configura- 
tions viz., 35 X 30, 30 x 26, 28 x 20 and 20 x 17, along z and r, respectively. In each 
case, the location of prediction has been kept nearly the same. A comparison of numer- 
ically predicted concentration profiles of different grid distribution is shown in fig. (2. 5). 
The values of the corresponding predicted mixing time for different grid system, as sum- 
marised in table (2.3), indicates that both 30 x 26 and 28 x 20 grid systems give almost 
same value of mixing time. Consequently, to be on the safer side, numerical predictions 
via a 30 X 26 grid system is considered to be essentially grid independent. 

2.4.3 Influence of incremental step size on predicted mixing 
time 

In earlier sections, the unsteady liquid mixing equation has been coupled with steady fluid 
flow and turbulence equations to predict the mixing time in a gas stirred bath. However, 
due to two dimensional nature of the transient mass conservation equation, an iterative 
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Table 2.3: Predicted variation of mixing time as a function of grid size distribution. 


SL. No. 

Grid Size 

mixing time 
(second) 

Location 

(r,z)cm 

1 

30 x26 

69sec 

(54.8,1.43) 

2 

35 x30 

77sec 

(53.0,3.63) 

3 

28 X 20 

69sec 

(54.4,1.53) 

4 

20 X 17 

152sec 

(54.1,2.22) 


d 

o 
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Figure 2.5: Variation in predicted mixing time as a function of grid size distributions ( 
At = 5sec, L = 0.8m, dt = 0.015m, gas flow rate = 6.67 x lO"** m^/s, and Cb = 0.0) 
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scheme has to be employed to obtain the normalized concentration profile numerically. 
The discretisation equation of the unsteady mass conservation equation includes time step 
(At) term (Eq. 2.24). Like grid size distribution, numerical predictions are also sensitive 
to time step. Hence, it is also important to establish solutions independent of the time 
step for further analysis. For this, normalized concentration profiles in a axisymmetric 
gas stirred cylindrical vessel (L=0.^0m and R=0.56m) have been predicted numerically 
by applying four different time steps viz., 3.5, 5, 10 and 20 seconds respectively. All 
otlior i)aramf!tcrs worn, koj)t fixed. Figure (2.G) represents the normalised concentration 
profiles for different time step values. The corresponding values of predicted 95% bulk 
mixing time for different time step (at the same location in all cases ) have been presented 
in Table (2.4). It is clear from Table (2.4) that both 3.5 and 5 seconds produce same 
values mixing time. Hence on the basis of above analysis, 5 second is considered as a 
reasonable value of time step (At) for numerical simulation. 

Table 2.4: Predicted variation of mixing time as a function of different time step at a 
particular location of a gas-stirred bath 


SL. No. 

Time step 
(second) 

Mixing time 
(second) 

Location 

(r,z)cm 

1 

3.5 

73sec 

(52.5,1.43) 

2 

5.0 

73sec 

(52.5,1.43) 

3 


93sec 

(52.5,1.43) 

4 

20.0 

139sec 

(52.5,1.43) 


2.4.4 Comparison with experimental results 

To assess the predictive capability of the mixing time model, normalized concentration 
profiles in gas stirred bath for different liquid depths viz., 0.5, 0.6, 0.8, 0.93 and 1.0 meter 
for a fixed vessel radius = 0.56m have been i)rodicted numerically and from these profiles, 
95% bulk mixing times were calculated for each configuration. Due precaution was taken 
so that the predicted mixing time truly reflects the bulk mixing time. Table (2.5) shows 
the value of the mixing time as function of liquid depth in the vessel. This indicates that 
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Figure 2.6: Variation in predicted mixing time as a function of four different incremental 
time steps (3.5, 5, 10 and 20 seconds), illustrates the specific value of At for time step 
independent numerical prediction ( grid Size = 30 x 26, L = 0.8m, R =0.56m, dj, = 
O.OlSra, gas flow rate = 6.67 x 10“'^ m^/s, C5 = 0.0) 
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mixing time decreases as depth of liquid increases. This is expected because the potential 
energy input rate of the system is directly proportional to the to the depth of liquid. In 
Fig. (2.7) the theoretically predicted mixing time has been plotted as a function of liquid 
depth and compared directly with experimental measurement of Mazumdar and Guthrie 
[5]. Reasonable agreement between prediction and measurement is readily apparent. 

Table 2.5; Predicted variation of mixing time as a function of liquid depth of a gas-stirred 
bath 


SL. No. 

Bath height 
(meter) 

Mixing time 
(second) 

Location 

(r,z)cm 

1 

0.5 

93 

(52.5,4.5)) 

2 

0.6 

80 

(52.5,5.3) 

3 

0.8 

71 

(52.5,4.3) 

4 

0.93 

48 

(52.5,4.98) 

5 

1.00 

43 

(52.5,5.35) 


2.4.5 Influence of empirical model parameters on predicted mix- 
ing time 


The interphase friction forces {F^andFr) and the turbulence production (Pj) due to bub- 
bles (Eq. (2.21) and Eq. (2.22)), are dependent on parameters such as bubble diameter 
(db), turbulence production constant (Cb) etc. The turbulence constant (Cj) takes a 
value between 0 and 1, and this is a function of bubble population, shape and size etc. It 
is obvious that the values of these parameters applied in the numerical solution scheme 
would affect the mixing time. While in earlier calculations the mean experimental value 
of bubble diameter was applied, additional calculations were carried out by considering 
different bubble diameters to assess the sensitivity of predicted results to 4- Thus a new 
bubble diameter has been calculated by applying following expression [29]: 


db 


0.35 



0.2 


(2.34) 




Time , 
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Liquid depth , m 

Figure 2.7: Theoretically predicted mixing times as a function of liquid depth for a 
gas-stirred vessel for( grid Size = 30 x 26, R = 0.56m, db = 0.015m, gas flow rate = 
6.67 X lO""* m^/s, Cj, = 0.0 ) and their comparison with experimental results. 
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Where; 

Q = Gas flow rate, 
g = Acceleration due to gravity. 

The value of dt, calculated from above equation is approximately 0.01 m. In the subsection 
2.4.4, theoretical prediction of mixing times were carried out for dj = 0.015 (experimental 
value). Table (2.6) represents the values of mixing time for five bath depths viz., 0.5, 0.6, 
0.8, 0.93 and l.Oin, deduced for a dt, value of 0.01m. This shows that, predicted mixing 
time for dt, =0.01m is more close to the observation of Mazumdar and Guthrie[5] than 
that of db = 0.015m. It is to be emphasized that the exact value of bubble diameter 
used in the numerical solution scheme is likely to influence the predicted results to some 
extent. 

Subsequent to this, influence of Cb values on predicted mixing time has been analysed. 
Mixing times were predicted numerically for two different value of Cj viz., 0.5 and 0.7 and 
compared with experimentally measured results[5]. The comparison shows (Fig. (2.9)) 
that predicted mixing time is sensitive to the previous values of the empirical constant 
Cb. 

Table 2.6: Predicted estimates of mixing time for theoretically deduced bubble diameter 
value of 10mm. 


SL. No. 

Bath height 
(meter) 

Mixing time 
(second) 

Location 

(r,z)cm 

1 

0.5 

100 

(52.5,4.5)) 

2 

0.6 

87 

(52.5,5.3) 

3 

0.8 

78 

(52.5,4.3) 

4 

0.93 

61 

(52.5,4.98) 

5 

1.00 

56 

(52.5,5.35) 
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Liquid depth , m 


Figure 2.8: Theoretically predicted mixing times as a function of liquid depth, for dif- 
ferent bubble diameters (0.015 m and 0.01m) and their comparison with experimental 
observations) 




Time , 


2.4 Results and Discussion 


35 



Liquid depth , m 


Figure 2.9: Theoretically measured mixing time as a function of liquid depth for different 
values of turbulence production constant (Cb = 0.5and0.7) and their comparison with 
experimental Jesuits. 
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2.4.6 The spatial dependency of mixing: the influence of scale 
factor (the size of the vessel) 

As pointed out previously, some investigators have reported on the variation of mixing 
times at different locations in the vessel, while others have reported contradictory ob- 
servations. An analysis of these investigations show that vessel size as well as degree 
of mixing were different in the studies. Consequently, mixing time was predicted in a 
smaller size vessel (L=0.24m and R=0.35m) to primarily assess the influence of ’scale 
factor’ on the rate of mixing at various locations. In Fig. (2.10) variation of mixing 
rates in three different regions of the vessel are shown. The nature of the curves for 
respective regions are essentially identical to those aJready shown for a bigger size vessel 
under identical operating gas flow rates. In Figure (2.11), locad mixing times for the three 
regions have been plotted as a function of degree of mixing (varied arbitrarily between 
90 to 99% )for the two different vessels. It is observed that 

(i) the magnitude of difference between the mixing times at various locations tend to di- 
minish as the degree of mixing is made more stringent. 

(ii) Sraaller the vessel size, smaller is the difierence between local mixing times for any 
degree of mixing. 

This indicates that i)osition dependency of mixing is not universal and is a function of 
vessel size and the criterion of mixing applied. 



2.4 Results and Discussion 



dng time for three different regions, A, B and C in a 
= 12 X 18, L = 0.24m, R = 0.35m, dh — 0.015m, gas 

Cft = 0.0) 
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Degree of mixing , % 


Figure 2.11: Effect of degree of mixing (viz., 90%, 92%, 95%, 97.5% and 99%) on local 
mixing times in different type of gas stirred vessels 



Chapter 3 


PREDICTION OF MASS 
TRANSFER RATES BETWEEN 
SOLID AND LIQUID IN GAS 
STIRRED BATH 

3.1 The General Method of Prediction of Solid-Liquid 
Mass Transfer Rates 

Numerous mass transfer correlations have been developed for forced convective mass 
tnuishM- situation onibodying throe key dimonsionloss mnnbors vi?:., Sh, Ro, Sc. These 
mass transfer correlations generally assume a functional form of the type: 

Sh = f{Re, Sc, geometry) 

Using, such correlations, mass transfer rates between solid-liquid can be predicted either 
from a purely experimental or a combined experimental and theoretical stand point. In 
case of mass transfer from a solid sphere or cylinder, the Sherwood number (Sh) is related 
to dissolution rate ) while Reynolds numbers is related to the flow cliaracteristics. 
Under the turbulent flow conditions, contribution to mass transfer arises from nominal 
Reynolds number as well as turbulent Reynolds number. Mass transfer rates in an 
axisyrnmetric gas stirred bath therefore can be predicted via any available correlation, 
using measured flow fields or can be directly obtained by monitoring the dissolution rate 
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experiinentally. 


3.2 Methodology of Developing a Mass Transfer Cor- 
relation 

Most of the mass transfer correlations have been developed by curve-fitting to experi- 
mental data. Various techniques can be adopted to develop mass transfer correlations, 
include viz., curve-fitting (linear or multiple regression) to a set of experimental measure- 
ments, curve-fitting to a set of data derived partly experimentation and partly through 
theoretical analysis, dimensional analysis, approximate integral technique, exact solution 
technique, etc. In addition to the above, modern artificial intelligence techniques such 
as, GAs and Artificial Neural nets (ANN) can also applied to a given same set of data to 
optimize the functional form of a given correlation. In this chapter possible application 
of GAs is examined from the view point of obtaining a mass transfer correlations with 
optimal pre-exponent and exponents. 


3.3 Possible Application of Genetic Algorithms for 
Correlation Development 

GAs are robust search and optimization techniques and are applicable to solve comj)lex, 
multimodal, discrete or discontinuous problems. The details of technique are provided 
elsewhere [30,31,32]. Flow diagram of (SGA) code used in present work is provided in 
Fig.3.1. 
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3.4 Derivation of Correlation from a Set of Data via 
GAs 


Experimental data of Singh and Mazumdar [33] and combined theoretical and experi- 
mental data of Mazumdar and coworkers [17] on dissolution of vertical solid cylinders in 
an aqueous model of gas stirred ladle has been considered for developing the correlations. 
Mazumdar [18], has shown that mass transfer rates from vertical cylinders in gas agitated 
reactors can be adequately described via a relationship of the type: 

Sh = 0.73Relf^Rel^^Sc°-^ (3. 1) 

In the preceding equation, 

Sh = Sherwood number. (^) 

Reioc,r = Nominal local Reynolds number. — ) 

Ret = Turbulent Reynolds number. (^) 

It is important to mention here that above correlation was obtained by multiple linear 
regression analysis of the measured data on dissolution. This equation for the experi- 
mental configuration of Mazumdar (Sc=500) can also be represented as: 

Sh = (3.2) 

In the most general form, eq.(3.2) can be written as: 

Sh = Re£, Re^^ (3.3) 

Where, A'l is the pre-exponent term and are the exponent terms. The purpose 

is to determine Ki , K 2 and K 3 through a fresh analysis (by the application of GA) of 
the reported data. The objective function for minimization is the sum of the residual 
squares, 
e.g., 


± ({Sh), - (if. (itef*).))' 


(3.4) 
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w^here, 

n is the total iiuiuber of data points. The set of data are presented in Table(3.2) and 
Table(3.3) respectively. 

The numerical value of various GA parameters used in the present investigation are 

• P, = 40 

• Per = 0.95 

• P^ = 0.01 

• Number of generations = 100 

• Total string length (1) =20, (First variable = 8, second variable=6 and third variable 
= ()) 

• Appropriate range of variables, is shown in Table(3.1), the range has been selected 
by considering the large number of correlations reported in standard texts and 
relevant literature. 

Table 3.1: Values of variable bounds applied in GAs correlation prediction 


SL.NO. 



variables 

lower value 

upper value 

1 




2 


0.2 

0.8 

3 

Kz 

0.25 

0.8 


3.4.1 Determination of Ki, K 2 and Kz on the basis of experi- 
mental data only 

Optimal solution is obtained considering experimental data only. Total number of data 
(n) is 32. This is shown in Table 3.2. After first generation the optimal mass transfer 
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correlation obtained by GAs is: 


Sh = 6.264 (3.5) 

Where, pre-exponent term = AT, = 6.264, and exponent terms are: 

K 2 = 0.23, K 3 = 0.35 

Hence, for Sc = 500, the above equation can be rewritten as: 

Sh = 0.80 (3.6) 


3.4.2 Determination of Ki, K 2 and Kz on the basis of theoretical 
and experimental data 

Optimal solution in this case is obtained by considering a set of combined experimental 
and theoretical data. Total number of data used (Table 3.3) is 56. The optimal mass 
transfer correlation by GA after first generation is 

Sh = 6.553 Re^ll Re^^'^ (3.7) 

Where the pre-exponent term Ki is found to be 6.553, and exponent terras are respec- 
tively K 2 = 0.22, and K 3 = 0.34. Hence, for Sc = 500, Equation (3.7) can also be written 
as: 


Sh = 0.84 Re^^ll Re°-^ (3.8) 

It is interesting to note from the above that regardless of the nature of data (whether 
experimental or theoretical) applied to GAs, correlation is practically identical results 
(see Eqs. (3.5) and (3.7)). This essentially indicates that the theoretical data which 
were obtained through the numerical solution of the turbulent Navier-Stokes equations, 
match well with experimental data. This implicitly indicates that mass transfer rates 
can be predicted from first principles with reasonable certainty. Since the GA based 
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Table 3.2: The set of experimental data of dimensionless numbers [33] considered for 
developing a mass transfer correlation by Genetic Algorithm. 


mm 

Sh 

Rei 

Rct 

1 

170.8 


93.9568 

2 

164.4 


95.0656 

3 

160.6 

429.05 

96.1072 

4 

11^ 

440.51 

126.0846 

5 


447.48 

128.0772 

6 

182.4 

453.99 

129.9402 

7 

160.28 

356.384 

128.871 

8 



130.4586 

9 

166.2 

365.12 

132.03 

10 

156.8 

529.587 

118.827 

11 


537.529 

120.609 

12 

161.4 

545.53 

122.40396 

13 

156.60 

626.19 

163.0524 

14 

153.16 

570.596 

148.5764 

15 

iQRai 

575.939 

149.9676 

16 

167.8 

689.068 

162.516 

HHBI 


673.651 

158.88 



680.435 

160.48 

19 

163.94 

466.554 

164.86 



449.63 

158.88 



453.705 

160.32 

22 

147.2 

754.3 

150.86 

23 

159.6 

792.6 

158.52 

24 

150.98 

800.36 

160.072 

25 

143.6 

660.346 

161.06 

26 

146.32 

666.004 

162.44 

mm 


944.809 


28 

154.10 

953.438 


■Bl 


685.608 

173.019 


138.90 

691.374 

174.4742 

31 

145.0 

808.05 

161.61 

32 

143.08 

815.5 

163.1 
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correlations are identical to the original correlation of Mazuindar, multiple regression is 
equally suitable for this class of problem. 


Table 3.3: Theoretically estimated data [17] of dimensionless numbers applied to develop 
a mass transfer correlation by Genetic Algorithm 



mam 



1 


328.0 

212.0 

2 

IBB 

M 


3 


11^ 


4 

140.0 

186.0 

130.0 

5 


396.0 

232.0 

6 

162.2 

617.0 

287.0 

7 

322.05 

1002.0 

984.0 

8 


137.0 

260.0 

9 

170.0 

297.0 

252.0 

10 

221.8 





1172.0 

723.0 

12 

204.9 


227.5 

13 

151.8 

596.0 

200.0 

14 

141.95 

819.0 

113.0 

15 

328.2 

2101.0 

586.0 


IBEI 

488.0 

270.0 


IBBi 






140.0 



2791.0 

630.0 

20 

186.05 

647.0 

356.0 

hieii 

IBBHI 

948.0 


mam 

mil 

763.0 

EBD 



HI 

EE 

24 

to 

^11 
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3.5 Prediclioiis of Moss Transfer Rates Tlirou^li 
Numerical Solution of Turbulent Navier-Stokes 
Equations 


The distribution of flow variables (u and v) and turbulence parameters («, e,ti etc.) in 
axisymmetric gas stirred cylindrical vessels (see Table (3.4)) have been predicted the- 
oretically by the solution of turbulent Navier-Stokes equations in conjunction with the 
K—e turbulence model so as to predict the solid-liquid mass transfer rates. The governing 
equations together with the set of boundary conditions and the numerical solution proce- 
dure have already been explained in Chapter 2. The velocity components and turbulence 
kinetic energy in two different size vessels (see Table 3.4) were predicted. Prom the pre- 
dicted distribution the nominal Reynolds number and the turbulent Reynolds numbers 
were estimated for those locations where the solid samples were kept submerged. To 

Table 3.4: Experimental conditions used for computation of mass transfer rates through 
numerical solution of of turbulent Navier-Stokes equation. 


Case 

Bath Depth 

Bath Radius 

Gas Flow Rate 

dfj 

Ct 


Vessel: 1 

0.21m 

0.15ra 

1 1/min 

O.Olm 

0.7 

10 X 8 

15 X 10 
20 X 16 

Vessel: 2 

0.42m 

0.238m 

8 l/rnin 

0.015m 

0.7 

24 X 20 
20 X 16 

15 X 10 


Table 3.5: Summary of thermophysical properties employed in the present work. 


Diffusivity of 


benzoic acid 

2.0 X 10-® 

in water m^/s 


Density of water, kg/rn^ 

1.0 X 10^ 

Viscosity of water, kg/m.s 

1.0 X 10-3 


I 


this end, the resultant mean velocity components, + were embodied in the ex- 
























3.5 Predictions of Mass Transfer Rates Tlurough Numerical Solution of Turbulent 
Navier-Stokes Equations 


48 


pression of nominal Reynolds number while the isotropic fluctuating velocity components 
calculated from u = were used in the expression of turbulent Reynolds number. 

The resultant velocity components together with the corresponding isotropic fluctuat- 
ing velocity components, thus obtained, are summarized in Table 3.6 and Table3.7 as a 
function of four different locations and three different grid systems (for the two vessels 
condition). These also include the values of nominal Reynolds number, Reioc,r 
and turbulent Reynolds number. Ret {^) for all configurations for three different grid 
systems. The numerical values of thermophysical properties employed to calculate above 
dimensionless numbers are presented in Table 3.5. 


Table 3.6: Predicted distribution of velocity fields and dimensionless numbers for axi- 
symmetric gas-stirred vessel ( L=0.42,R= 0.238 and d=10mm ) 


grid 

size 

loca 

tion 

z 

III 

r 

m 

y/u^ + V^ 

u 

Re 

loc 

Ret 

Sh 

exp 

Sh 

pred 


1 

0.228 

0.104 

0.0273 

0.0198 

273 

198 

112.7 

127.6 

24 X 20 

2 

0.06 

0.046 

0.0189 

0.029 

189 

290 

131.65 

131.5 


3 

0.383 

0.047 

0.1147 

0.1261 

1147 

1261 

283.9 

330.4 


4 

0.357 

0.106 

0.0466 

0.0672 

466 

672 

140.0 

215.6 


1 

0.228 

o.ioP 

0.032 

0.0533 

320 

533 

112.7 

182.3 

20 X 16 

2 

0.06 

0.046 

0.0243 

0.033 

243 

330 

131.65 

145.9 


3 

0.383 

0.047 

0.1145 

0.1634 

1145 

1634 

283.9 

358.8 


4 

0.357 

0.106 

0.0517 

0.0972 

517 

972 

140.0 

249.0 


1 

0.228 

0.104 

0.0313 

0.1065 

313 

1065 

112.7 

226.5 

15 X 10 

2 

0.06 

0.046 

0.0424 

0.0875 

424 

875 

131.65' 

229.2 


3 

0.383 

0.047 

0.065 

0.1897 

650 

1897 

283.9 

326.7 


4 

0.357 

0.106 

0.034 

0.1361 

340 

1361 

140.0 

249.8 


3.5.1 Estimation of grid size independent solutions 

I 

Theoretically predicted flow variables and turbulence parameters in gas stirred vessel are 
expected to be sensitive to the grid size distribution. It is important to assess the influence 
of grid size on predicted results. Sherwood numbers have been calculated by using the 
correlation of Mazumdar et. al. [18] {Sh = for all the conditions and 



3.5 Predictions of Mass Transfer Rates Through Numerical Solution of Turbulent 
Navier-Stokes Equations 


49 


Table 3.7: Predicted distribution of velocity fields and dimensionless numbers for ax- 
isyrametric gas stirred vessel ( L=0.21 and R= 0.15 ) 


grid 

size 

loca 

tion 


Bi 

mm 


Re 

loc 

Ret 

Sh 

exp 

Sh 

pred 


1 


0.095 

0.00884 


150 


145 

157.7 

.10 X 8 

2 


0.055 



458.37 

1173.43 

143 

256.7 


3 

0.04 

0.035 



273.76 

717.27 

133 

192.8 


4 

0.14 

0.045 

0.0285 

0.08 

424.36 



253.0 


1 




0.029 

183.3 

492.18 

145 

154.6 

15 X 10 

2 


0.055 


0.058 

395.15 

916.75 

143 

228.6 


' 3 


0.035 


0.0373 

349.5 

600.68 

133 

193.6 


4 

0.14 

0.045 

0.036 

0.0672 

536.04 



253.7 


1 


0.095 

0.0078 


132.38 

147.65 

iim 

96.95 

20x 16 

2 


0.055 

0.0174 



245.0 


136.86 


3 


0.035 

0.0115 


185.2 

95.0 

133 

91.56 


4 

0.14 

0.045 

0.0297 


442.2 

376.7 

130 

176.86 


compared with corresponding experimental values of Sherwood number. Tables 3.6 and 
3.7 provide values of experimental and predicted Sherwood number. It is seen that for the 
finest grid system only the experimental and predicted values are comparable. Hence, the 
results deduced via the finest grid systems (20 x 16 for first vessel (L=0.21,R=0.15) and 
24 X 20 for second vessel (L=0.42,R=0.238)) have been considered for further analysis. 

3.5.2 Comparison of numerical prediction with GAs-based cor- 
relation 

Table 3.8 represents the numerically predicted dimensionless numbers viz., Sh, Reioc,r, 
and Ret for the optimal grid systems which are essentially reproduced from the previ- 
ous tables, lb assess the adequacy of Eq. (3.6) (mass transfer correlations developed 
by GAs using experimental data) with reference to the set of numerically predicted re- 
sults , a log-log plot between the quantities Sh/Re^^^ and Re^^ is been shown in Fig. 
(3.2). The straight line corresponding to the GAs correlation (Eq. (3.6)) is shown 
together with the numerically predicted results. It is clear that GAs based correlation 
{Sh = 6.264Re|’^fr^6P®) agrees well with the numerically predicted points. Finally, 
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Table 3.8: Values of dimensionless groups obtained via numerical solution of turbulent 
Navier-Stokes equation for an optimal grid configuration 


SI. No. 

R^LoCfT 

Ret 

Sh 



Sn 

„Loc - 

1 

273 

198 

112.7 

32.8 

31.0 

2 

189 

290 

131.65 

41.55 

39.43 

3 

1147 

1261 

283.9 

60.265 

56.165 

4 

466 

672 

140.0 

36.23 

34.07 

5 

132.38 

147.65 

145.0 

49.5 

47.13 

6 

275 

245 

143.0 

41.56 

39.3 

7 

185 

95 

133.0 

42.17 

40.02 

8 

442.2 

376.7 

130.0 

34.03 

32.02 


t/3 



CP 


on 

o 


4 

3 

2 

1 

0 


1 1 1 1 1 1 r 

Present work ( L=0. 42m , R=0.238in) ^ 
Present work ( L=0.21m , R=0.15m ) + 
Present work (GAs correlation ) 



I I I I I 1 1 1 1 

0 1 2 3 4 


Jog Rej 


Figure 3.2: Logarithmic variation of and Ret illustrating the adequacy of GAs mass 
transfer correlation with reference to tlie numerically predicted mass transfer rates. 

the adequacy of Eq. (3.8) (mass transfer correlation developed by GAs using both ex- 
perimental and theoretical data Table 3.3) has been assessed in a similar way to that 
mentioned above. A log-log plot between the quantities Sh/Re^g^ f. and Ret has been 
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shown in Fig. (3.3). This also indicates close agreement between GAs based correlation 
and numerically predicted results. 



log Re^ 


Figure 3.3: Logarithmic variation of and Ret illustrating the compatibility of GAs 

Loc 

mass transfer correlation with reference to the numerically predicted mass transfer rates. 





Chapter 4 


CONCLUSIONS 


The following major conclusions can be drawn from the present investigation regarding 
mixing time prediction and mass transfer correlation. 


4.1 Mixing time prediction 

• Mixing rates in gas stirred melts have been theoretically predicted and compared 
with available experimental results. It is observed that mixing phenomena in such 
systems can be predicted from the first principles with reasonable accuracy. 

• The region close to the bottom of ladle, where flow is relatively weak, exhibited the 
slowest rate of mixing in the system. This agrees with earlier work [5]. 

• Numerical results shows that as the liquid depth in the vessel increases the mixing 
time decreases. The same phenomena is observed in experiments. 

• Mixing time depends on the value of empirical constants (dj, Cb etc.) applied. 
Hence, for realistic prediction sufficiently accurate information on these are re- 
quired. 

• The spatial variation of mixing is a function of the scale of the vessel and the mixing 
criteria adopted. 
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4.2 Mass transfer correlation 

• The adequacy of a previously reported mass transfer correlation has been re- 
examined by the application of Genetic Algorithms. A large number of data 
reported in the literature on dissolution of vertically submerged solids has been 
analysed, and an optimal set of values for the pre-exponent and exponent terms 
derived via GA. The resultant correlation is practically identical to the one deter- 
mined earlier by Mazumdar et al. by using linear regression. 

• Predicted mass transfer rates via the GA based correlation was found to agree well 
with numerical solution of turbulent Navier-Stokes equation. 



Chapter 5 


RECOMMENDATIONS FOR 
FUTURE WORK 


The scope of the present work can be further expanded by including the following: 

• Mixing time model developed in the present work may be applied to investigate 
mixing in other configurations of gjus stirred bath viz., CAS etc. 

• Mass transfer correlation may be assessed against industrial data on dissolution of 
ferro-alloys etc. 

• Application of combined GAs-ANN method remains to be explored to simulate 
mass transfer from theoretical and industrial data. 
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